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(5)  INTRODUCTION 

Subject:  This  study  concerns  the  development  of  MRI  arterial  spin  tagging  to  non-invasively 
measure  breast  tissue  perfusion.  The  project  will  test  whether  the  tagging  technique  is  a  low  cost 
and  more  reliable  alternative  to  dynamic  first-pass  contrast-enhanced  imaging,  giving 
possibilities  for  use  in  specific  clinical  settings,  as  well  as  in  screening.  MR  dynamic  first-pass 
contrast-enhanced  imaging  has  shown  that  malignant  and  benign  breast  lesions  can  be 
distinguished.  However,  it  may  have  limited  importance  in  clinical  breast  diagnosis  due  to 
significant  false-negatives  and  false-positives.  Also,  dynamic  contrast  enhancement  cannot  be 
used  in  screening,  even  within  defined  high-risk  sub-populations  such  as  patients  with  dense 
breasts,  due  to  the  high  cost  of  contrast  material  and  administration.  The  arterial  spin  tagging 
technique  was  developed  to  measure  tissue  perfusion  parameters  without  the  use  of  contrast,  and 
has  been  successfully  demonstrated  in  brain  and  kidney.  The  spin  tagging  technique  may  provide 
more  sensitive  and  specific  assessment  of  these  parameters,  at  higher  spatial  resolutions,  by 
virtue  of  its  ability  to  acquire  repeated  measurement  of  flow  sensitive  data  to  improve  signal  to 
noise  ratio  (SNR).  Dynamic  first-pass  contrast-enhanced  imaging  has  a  limited  data  acquisition 
time  and  thus  a  limited  resolution  and  SNR. 

Purpose:  The  purpose  (long-term  goal)  of  the  project  is  to  develop  MRI  arterial  spin  tagging  as  a 
low  cost  and  more  reliable  alternative  to  first-pass  contrast-enhanced  imaging.  We  will  test  the 
hypothesis  that  arterial  spin  tagging  provides  accurate  and  precise  discrimination  between 
normal  tissue,  benign  and  malignant  lesions,  when  differences  in  perfusion  and  Ti  exist.  Lesions 
will  have  been  previously  detected  by  clinically  accepted  diagnostic  imaging  procedures,  and  by 
biopsy.  Statistical  analysis  will  be  performed  to  access  the  correspondence  between  arterial  spin¬ 
tagging  and  biopsy,  and  to  establish  the  relative  value  of  spin  tagging  compared  to  first-pass 
contrast-enhanced  MRI.  We  hope  to  establish  that,  relative  to  using  first-pass  contrast-enhanced 
imaging,  false  positives  and  negatives  are  reduced  using  arterial  spin  tagging  by  virtue  of 
increased  image  signal-to-noise  ratio  (SNR),  higher  spatial  resolution,  and  the  unique  ability  to 
obtain  estimates  of  macromolecular  bound  fluid  fractions. 

Scope:  The  effort  on  the  project  is  mainly  limited  to  the  technical  aspects  of  development  of  a 
new  methodology.  However,  the  project  also  includes  a  rigorous  performance  comparison  with 
the  current  gold-standard  methodology.  The  specific  aims  are  to  (1)  refine  arterial  spin  tagging 
pulse  sequences  and  imaging  protocols,  (2)  develop  automated  data  analysis  software  for 
measurement  of  breast  tissue  parameters,  and  (3)  compare  the  technique  to  first-pass  contrast- 
enhanced  MRI  and  biopsy.  The  technique  will  be  evaluated  in  sixty  patients,  with  roughly  equal 
numbers  of  benign  and  malignant  lesions.  Multi-component  differential  equation  models  of 
magnetization  distribution  within  breast  tissue  will  be  used  to  estimate  pension  and  Tu  from 
which  regions  suspicious  for  malignant  disease  will  be  identified.  While  sub-population 
screening  is  a  critical  potential  application,  it  is  premature  to  propose  a  study  to  evaluate  the 
technique  for  specific  sub-population  screening,  before  it’s  capability  for  detection  and 
discrimination  of  previously  detected  lesions  is  established. 

Background:  Studies  have  shown  that  dynamic  first-pass  contrast-enhanced  magnetic  MRI 
discriminates  malignant  from  many  types  of  benign  breast  disease.  This  technique  detects 
malignant  disease  by  measuring  tissue  fluid  distribution  parameters,  which  are  derived  from  the 
rate  of  rise  of  MRI  signal  intensity.  The  rate  of  decrease  of  signal  intensity,  after  the  peak  signal 
is  reached,  is  also  believed  to  be  related  to  malignancy,  based  on  the  idea  that  the  vasculature  of 
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malignant  tissues  possesses  a  high  degree  of  arteriovenous  shunting  causing  more  rapid 
elimination  of  contrast  agent.  Unfortunately,  this  technique  may  have  limited  importance  due  to 
significant  false-negatives  and  false-positives.  Further  development  and  refinement  is  required, 
but  major  improvements  are  unlikely  given  the  limited  data  acquisition  time  inherent  to  first-pass 
techniques.  In  particular,  the  first-pass  techniques  cannot  provide,  with  high  temporal  resolution, 
high  spatial  resolution  required  to  reveal  morphological  details  that  can  reduce  false  positives. 
Arterial  spin  tagging  has  been  developed  for  measuring  tissue  parameters  without  the  use  of 
contrast,  and  has  been  successfully  demonstrated  in  brain  and  kidney.  Spin  tagging  measures 
absolute  blood  perfusion  (cc/min/g  tissue),  longitudinal  relaxation  time  (T\)  of  the  tissue  water, 
and  possibly  macromolecular  bound  fluid  fractions.  Spin  tagging  can  provide  more  sensitive  and 
specific  assessment  of  tissue  parameters,  at  higher  spatial  resolution,  using  the  ability  to  build  up 
signal  to  noise  over  reasonable  scan  times  to  investigate  tissue  parameters.  It  is  not  limited  to 
data  collection  during  a  first  pass. 

Previous  studies  have  confirmed  that  first-pass  enhancement  of  carcinomas  by  Gd-DTPA  are 
often  (but  not  always)  more  intense  and  faster  compared  to  normal  parenchyma  or  benign  breast 
disease  [1,2,3,4,5,6,7,8].  Enhancement  is  attributed  to  tumor-induced  angiogenesis  resulting  in 
increased  vascularity  within  the  malignant  tumor  [9,10,1 1]  which  often  also  develops 
anastomoses  and  shunts.  Due  to  greater  tissue  uptake,  carcinomas  enhance  maximally  within  1.5 
to  2  minutes  of  contrast  injection  and  then,  due  to  accumulation  in  extracellular  fluid,  tend  to 
maintain  the  enhancement  for  several  minutes.  Kaiser  [12]  postulated  that  false  positive  benign 
lesions  may  be  distinguished  by  the  behavior  of  signal  enhancement  after  the  peak  occurs,  i.e., 
during  the  elimination  phase  [13].  Malignant  lesions  have  greater  degrees  of  AJV  shunting,  and 
hence  more  rapid  elimination  of  contrast  and  more  rapid  return  to  non-enhanced  signal  levels. 
Due  to  the  importance  of  the  time  course  of  signal  enhancement  in  these  criteria, 
pharmacokinetic  analysis  [8,14,15]  has  been  clinically  important.  Motion  correction  [16]  and 
image  subtraction  [17]  have  been  important  for  reliable  temporal  processing.  However,  the 
importance  of  morphology  is  regarded  by  many  researchers  as  equally  or  more  important  than 
temporal  profiles,  and  some  have  emphasized  high  spatial  resolution,  not  temporal  resolution,  in 
discrimination  of  tumors  [18,19  page  595].  Overall,  inconsistent  reports  and  general  overlap  of 
benign  and  malignant  temporal  profiles  precludes  the  role  of  MRI  as  an  “imaging  biopsy” 
[19,20]. 

Techniques  for  non-invasive,  spatially  resolved  MRI  perfusion  measurement  have  been 
developed  based  on  magnetic  labeling  of  arterial  water  [21,22]  (i.e.  arterial  spin-tagging),  and 
tested  in  both  animals  and  humans  [23,24,25,26].  Mathematical  models  of  the  magnetization 
distribution  are  used  to  relate  changing  signal  intensity  to  tissue  perfusion  (flow  of  magnetization 
into  tissue),  magnetization  of  free  fluid,  magnetization  of  fluid  bound  to  macromolecules,  and 
fluid  exchange  rates  between  these  compartments  [23,27].  Tagged  spins  will  not  accumulate 
significantly  in  tissues  as  does  contrast  agent,  and  therefore  do  not  obtain  precisely  the  same 
measurements  as  contrast-enhanced  scans.  However,  arterial  spin  tagging  will  measure 
macromolecular  bound  fluid  fraction,  that  is  affected  by  extracellular  fluid  fraction,  and  can  be 
adapted  to  measurement  of  fluid  elimination. 
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(6)  BODY 


Overview 

This  section  details  the  experimental  methods,  results  and  discussion  in  relation  to  the  Statement 
of  Work  outlined  in  the  proposal.  The  first  subsection  provides  brief  paragraphs  of  the 
accomplishments  in  relation  to  the  Statement  of  Work.  Detailed  discussion  of  these 
accomplishments  is  provided  in  the  second  subsection.  There  are  two  main  contributions  of  the 
work.  The  first  is  the  discovery  of  a  multitude  of  unrecognized  sources  of  T\  estimate  errors,  the 
recognition  of  the  profound  sensitivity  of  perfusion  measurements  to  T\  estimate  errors,  and  the 
development  of  robust  pulse  sequences  and  analysis  methods  to  derive  accurate  T\  estimates  with 
inherently  noisy  and  confounded  signal  data.  The  second  is  the  development  of  the  software 
program,  BrView,  to  allow  rapid  review,  quantitative  analysis,  and  assessment  of  the  multitude 
of  different  breast  images  and  timeseries  data  that  is  obtained  in  each  patient  study. 

Summary  of  accomplishments  in  relation  to  Statement  of  Work 
Technical  objective  1 

Task  1:  Months  1-6:  Implementation  and  testing  of  magnetization  transfer  pulses  for  both  arterial 
tagging  and  first-pass  contrast  enhanced  sequences. 

Magnetization  transfer  pulses  have  not  been  implemented  yet.  Mainly,  this  delay  was  due  to  the 
fact  that  the  MRI  system  at  UC  Davis  Medical  Center  is  scheduled  to  be  upgraded  to  the  new 
Horizon  LX  system  soon,  and  we  opted  to  defer  this  pulse  sequence  development  until  we  were 
using  this  new  LX  software  platform.  There  have  been  multiple  delays,  but  this  upgrade  is 
scheduled  to  occur  by  January  1999.  The  existing  pulse  sequence  will  be  converted  from  the  5.4 
OS  platform  to  the  LX2  platform  (a  substantial  change),  and  the  magnetization  transfer  pulses 
will  be  implemented  also. 

Because  of  concerns  regarding  the  ultimate  utility  of  the  magnetization  transfer  technique,  we 
have  decided  that  implementation  and  testing  of  arterial  spin-tagging  sequence  based  on  echo 
planar  imaging  (EPI)  data  acquisition  is  a  higher  priority.  EPI  based  spin-tagging  technique  can 
theoretically  improve  the  accuracy  of  the  Ti  estimation,  which  we  have  studied  and  worked  with 
extensively  (see  below).  The  major  problem  with  magnetization  transfer  is  that  it  represents  a 
perturbation  on  an  already  small  signal.  Therefore,  whether  we  will  observe  a  magnetization 
transfer  effect  is  questionable.  Nevertheless,  it  will  be  implemented  as  part  of  the  project. 

Task  2:  Months  1-9:  Implementation  and  testing  of  interleaved  high-resolution  imaging 
technique,  for  both  arterial  tagging  and  first-pass  contrast  enhanced  sequences. 

An  interleaved  arterial  spin-tagging  technique,  that  allows  high  resolution  (e.g.  256  x  256  or 
more)  has  been  implemented  and  tested.  We  also  implemented  a  rectangular  data  acquisition 
technique,  which  acquires  data  sets  that  have  different  numbers  of  points  along  the  frequency 
and  phase  encode  direction.  The  256  x  240  interleaved  acquisition  provides  the  best  spatial 
versus  temporal  resolution  tradeoff,  and  has  been  used  for  all  of  our  recent  studies. 

Task  3:  Months  1-9:  Analysis  of  spin  tagging  sequence  to  understand  causes  of  existing  baseline 
offsets,  effects  of  inversion  slice  transition  profiles,  and  effects  of  RF  flip  angle  profile  on  the 
measurements,  with  implementation  and  testing  pulse  sequence  modifications  to  minimize  these 
imperfections. 
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Matching  of  inversion  and  excitation  slice  profiles  was  greatly  improved  by  using  customized 
RF  pulses  designed  with  the  Shinnar-Le  Roux  (SLR)  algorithm.  The  new  data  processing 
technique  (based  on  a  semi-log  linear  regression  of  inversion  time  {Tl)  dependent  signal) 
significantly  reduced  the  previously  reported  problem  in  measuring  the  longitudinal  relaxation 
time  (Tl)  caused  by  baseline  offsets,  and  by  the  uncertainty  in  the  effective  inversion  time.  The 
optimal  RF  flip  angle  for  spin  excitation  during  spatial  encoding  was  found  to  be  10  degrees, 
based  on  experiments  over  a  range  of  flip  angles. 

Technical  objective  2 

Task  4:  Months  :^-15:  Software  for  automatic  estimation  and  error  analysis  of  perfusion,  tissue 
water  longitudinal  relaxation  time,  and  extracellular  fluid  volume  fraction  from  mathematical 
models  and  user-defined  ROIs  from  spin  tagging  timeseries.  Implementation  of  pharmacokinetic 
model  calculations  and  error  analysis  for  first-pass  contrast  enhanced  imaging. 

A  software  program  (BrView)  has  been  written  in  C,  X  Window  System,  and  Motif,  and 
implemented  on  an  SGI  02  computer  (paid  for  by  grant  funds)  for  analyzing  and  visualizing  the 
MR  images  and  timeseries.  Much  effort  has  been  expended  on  developing  the  capability  to  easily 
review  all  of  the  images  taken  in  each  study,  and  cross-reference  pixel  locations.  Much  more 
effort  than  anticipated  has  been  expended  in  the  development  of  a  robust  technique  for 
identifying  so-called  “suspicious  regions”  based  on  the  T\  and  perfusion  measurements.  The 
measurement  of  Ti,  upon  which  the  perfusion  measurement  entirely  depends,  is  based  on  a  semi¬ 
log  linear  regression  technique  developed  by  the  investigators.  The  Ty,  perfusion  (f/X)  and 
standard  errors  of  these  quantities  are  estimated  automatically  using  this  robust  technique.  The 
so-called  feature  images  display  these  quantities  in  color  and  grayscale.  These  quantities  are  used 
to  calculate  a  “suspicion  index”  for  each  pixel,  and  thereby  identify  regions  of  breast  tissue 
suspicious  for  malignancy. 

The  implementation  of  the  pharmacokinetic  model  calculations  and  error  analysis  for  first-pass 
contrast  enhanced  imaging  has  not  been  completed.  Thus  far,  the  software  allows  the  user  to 
click  on  the  reference  image  of  the  dynamic  study,  to  display  the  time  profile  of  the  signal  at  that 
pixel.  The  implementation  of  the  pharmacokinetic  model  is  currently  a  major  focus  of  the  current 
effort. 

Task  5:  Months  3-15:  Software  for  registration  (including  implementation  and  testing  of  motion 
correction  and  physiological  noise  reduction  algorithms)  and  overlay  of  images  from  high 
resolution  Ti  weighted,  spin-tagging,  and  contrast  enhanced  studies. 

Using  the  BrView  software  program,  image  pixels  containing  high  values  of  perfusion  (f/X),  and 
moderate  to  high  values  of  Ti,  and  low  standard  errors,  are  identified  as  suspicious,  and  assigned 
a  “suspicion  index”  based  on  finding  similar  abnormal  values  at  spatially  adjacent  pixels. 
Suspicious  pixels  can  be  overlaid  on  the  high-resolution  T\  and  T2  clinical  images  by  color 
mapping,  and  simultaneously  presented  with  the  first-pass  contrast  enhanced  timeseries  at  these 
pixels.  This  provides  the  user  with  a  comprehensive  anatomical  and  functional  view  of  the 
suspicious  regions,  and  facilitates  making  a  decision  regarding  the  malignant  nature  of  the  lesion. 
Finally,  a  motion  artifact  estimation  and  correction  algorithm  has  been  developed  and 
implemented  in  a  separate  software  program. 
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Technical  objective  3 

Task  6:  Months  9-24:  60  patients  with  malignant  and  benign  breast  lesions  will  be  imaged  usin^ 
Ti  weighted,  first-pass  contrast  enhanced,  and  arterial  spin  tagging  MRI  pulse  sequences. 


We  have  recruited  only  two  subject  studies  for  both  first-pass  contrast  enhanced  and  arterial  spin 
tagging  MRI  pulse  sequences.  Now  that  the  software  development  is  substantially  completed,  we 
are  trying  to  improve  the  rate  at  which  patients  are  recruited  into  the  study.  It  is  likely  that  the 
patient  studies  will  extend  into  the  third  year  of  the  proposal.  We  are  currently  discussing  with 
physicians  at  a  local  Kaiser  Hospital  for  recruitment  of  additional  patients.  Appropriate 
documentation  will  be  submitted  for  approvals  upon  reaching  a  firm  agreement. 


Detailed  description  of  the  accomplishments 


Background  to  theory  of  arterial  spin  tas 


This  background  material  is  taken  from  the  literature,  and  represents  the  starting  point  for  the 
project.  In  the  so-called  “on-slice”  spin  tagging  technique  that  we  implemented,  spins  in  a  slab  of 
tissue  are  inverted  first  using  a  RF  pulse.  The  thickness  of  the  slab  can  be  controlled  by  the 
strength  of  the  slab  selective  gradient  (normally  a  Z  gradient).  Then  the  levels  of  recovery  to 


equilibrium  are  measured  at  different  times  (Figure  1).  The  rate  of  recovery  (i.e.  the  effective  Ti) 
can  be  estimated  from  the  recovery  curve.  In  the  research  protocol,  3mm  slices  are  imaged 
Magnetization  recovery  after  time  TI 


TI 


Figure  1.  After  the  Magnetization  is  inverted,  it  Figure  2.  A  voxel  within  the  imaged  slice  is 

recovers  towards  thermal  equilibrium  position  identified.  Arterioles  and  capillaries  pass  through  the 

according  to  characteristic  relaxation  time  Ti.  The  voxel  resulting  in  net  blood  flow  (indicated  by  long 
signal  intensity  along  a  typical  magnetization  recovery  arrows)  and  replacement  of  tissue  water.  Blo^  water 
curve  is  shown.  The  image  displays  the  magnitude  of  (indicated  by  short  arrows)  exchanges  readily  with 
the  magnetization.  The  rate  of  recovery  can  be  the  extracellular  fluid, 

determined  from  images  acquired  at  different  times  of 

the  recovery  curve  represented  by  TI.  TI  is  a  point  in  sequentially  (one  at  a  time).  The  goal  is  to 

time  at  which  an  zero  k-soace  line  is  acauired.  measure  capillary  blood  perfusion,  which  is 

defined  as  the  net  amount  of  arterial  blood 

per  second  entering  a  gram  of  breast  tissue.  This  definition  can  be  extended  to  a  per  unit  volume 
(voxel)  measurement,  i.e.,  the  net  amount  of  arterial  blood  per  second  entering  a  voxel  of  breast 
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tissue.  Therefore,  perfusion  represents  the  directed  microscopic  flow  within  a  voxel.  In  breast 
tissue,  blood  vessels,  which  are  basically  composed  of  capillaries,  are  very  small.  Blood  flow  is 
very  slow  in  these  vessels.  Water  molecules  can  move  easily  between  the  vessels  and  their 
surroundings  (Figure  2).  A  two-component  model  is  used  to  approximate  the  activity  of  the 
water  molecules.  In  this  model,  water  molecules  in  blood  vessels  entering  a  slice  of  breast  tissue 
mix  thoroughly  with  water  molecules  outside  the  vessels  before  they  leave  the  slice.  For  this 
two-component  model,  the  Bloch  Equation  for  the  breast  tissue  water  magnetization  at  the  slice 
of  tissue  investigated  can  be  written  as  [23] 


dM{t) 

dt 


Mq  -Mit) 

T, 


(1.1) 


where,  M{t)  =  the  longitudinal  water  magnetization  of  a  voxel  of  breast  tissue  at  time  t,  Mq  =  the 
value  of  M(0  under  fully  relaxed  conditions,  Ti  =  the  spin-lattice  relaxation  time  constant  of  the 
breast  tissue,  Mm(t)  =  the  longitudinal  magnetization  of  macromolecules  in  a  voxel  of  tissue,  kfar, 
krev  =  the  magnetization  transfer  rate  constants  from  water  to  macromolecules  and  vice  versa,  Ma 
=  the  longitudinal  water  magnetization  per  ml  of  arterial  blood  entering  a  voxel  from  outside  of 
the  slice  investigated, /=  perfusion  in  ml  of  arterial  blood  per  second  in  a  voxel,  and  X  =  the  ratio 
of  water  content  between  breast  tissue  overall  and  the  arterial  blood  within  it. 

To  simplify  Eq.(l.l),  we  assume  that  the  magnetization  transfer  contributions  are  negligible, 
specifically. 


\(Kev-kfor)M„(t)\ 


« 


Mq  -Mjt) 


T, 


+  -f 


Mit) 


Then  Eq.  (1.1)  can  be  simplified  to 

dt  r,  A  J 


(1.2) 


Eq.  (1.1)  will  be  utilized  without  this  approximation  when  magnetization  transfer  is 
implemented.  Parameters  Ti  mdf/X  in  Eq.  (1.2)  can  be  estimated  from  two  spin  tagging 
conditions.  In  the  so-called  “Selective  Tagging  Condition”,  all  spins  within  the  slice  of  breast 
tissue  being  imaged  are  inverted,  and  the  signal  from  the  recovering  magnetization  is  acquired  at 
a  specified  later  time  (defined  as  the  inversion  time  TI).  Images  are  acquired  at  several  TI  values 
in  order  to  generate  a  magnetization  recovery  curve  at  each  pixel.  Since  the  arterial  spins 
entering  the  slice  are  not  influenced  by  these  selective  inversion  pulses,  they  enter  the  slice  with 
full  longitudinal  magnetization  (i.e.  they  are  in  thermal  equilibrium).  Then, 

Ma(t)=Mf 


where  Ma  =  the  value  of  Mait)  at  thermal  equilibrium.  In  the  two-component  model,  we  assume 
that  the  arterial  blood  spins  entering  the  slice  mix  thoroughly  with  spins  in  the  surrounding  breast 
tissue  before  leaving  the  slice.  Thus, 
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and  Eq.  (1.2)  becomes 


(1.3) 


A  MR  system  measures  the  apparent  Tu  here  represented  as  T\^,  as  if  the  all  spins  are  stationary 
(not  flowing): 

dM{t)  _  Mo-M(0 
dt  ■  7;, 


where, 


(1.4) 


and  the  subscript  s  denotes  that  the  parameter  was  measured  using  the  selective  RF  pulse. 

Images  are  also  acquired  in  a  so-called  “Non-selective  Tagging  Condition”.  Here,  spins  are 
inverted  everywhere  in  the  sensitive  volume  of  the  body  coil,  but  signals  are  acquired  only  from 
the  slice  being  investigated,  as  in  the  Selective  Tagging  Condition.  The  arterial  spins  flowing 
into  the  slice  thus  experience  they  same  inversion  pulses  as  those  within  the  slice,  and  the 
magnetization  of  entering  arterial  blood  can  be  written 

/  M.{t)  =  /  ^ 


Eq.  (1.2)  reduces  to  the  simple  expression 


dM{t)  _ 
dt  ~  r, 


where  Ti  in  this  case  represents  the  true  T\  of  the  tissue  without  the  confounding  influence  of 
blood  flow.  This  Ti  and  corresponding  rate  b„  is  denoted 


and 


rin  =  7’i 


(1.5) 

(1.6) 


where  the  subscript  n  is  used  to  distinguish  these  parameters  from  those  found  in  Selective 
Tagging  Condition.  Based  on  the  b,  and  b„  values  obtained  from  the  two  tagging  conditions,  we 
can  determine  the  perfusion  parameter//A.: 


/ 

X 


=  b.  -K 


(1.7) 
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This //A,  is  the  rate  of  magnetization  replacement  in  a  voxel.  The  importance  off  /X  is  that, 
compared  to  normal  breast  tissue  and  benign  breast  lesions,  a  higher//A,  is  expected  in  malignant 
tumors,  due  to  higher  blood  flow  rates  into  the  malignant  tumor. 

Contributions  to  pulse  sequence  design 

The  University  of  California  Davis  Medical  Center  operates  a  Signa  Advantage  1.5  T  GE  MR 
system  (Milwaukee,  Wisconsin)  (with  1.0  g/cm  peak  gradient  strength,  500  microsecond  rise 
time,  and  version  5.4  OS)  for  MR  clinical  and  basic  science  research.  An  arterial  spin  tagging 
pulse  sequence  has  been  developed  for  this  MR  system  based  on  the  original  source  code  of  the 
commercial  GE  2D  fast  Spoiled  Gradient  Recalled  Echo  (SPGR)  pulse  sequence.  In  this 
sequence,  after  the  inversion  RF  pulse,  raw  data  is  acquired  for  an  entire  2D  image  using  the  fast 
SPGR  k-space  acquisition  (Figure  3).  Several  innovations  have  been  implemented  in  this  custom 
pulse  sequence:  1.  The  original  2D  fast  SPGR  data  acquisition  scheme  has  been  modified  to 
support  a  symmetric  centric  k  space  data  acquisition.  Data  acquisition  starts  at  the  zero  ky  line  of 
the  k  space,  and  alternates  below  and  above  this  ky  line,  end  at  the  highest  frequency  ky  lines. 
This  allows  the  low  spatial  frequency  components,  which  have  the  dominant  effect  on  image 
contrast,  to  be  acquired  at  the  beginning.  The  effective  time  of  inversion,  defined  by  the  time 
from  spin  inversion  to  the  acquisition  of  the  center  k  space  line,  can  be  adjusted  most  flexibly 
with  this  acquisition.  2.  Implementation  of  an  optimized  RF  inversion  pulse  and  an  optimized  RF 
acquisition  pulse,  both  derived  fi-om  the  Shinnar  Le-Roux  (SLR)  algorithm  [28].  These 
optimized  RF  pulses  allow  slice  profiles  with  sharp  edges  so  that,  in  the  selective  condition,  the 
alignment  between  the  inversion  slice  and  the  image  slice  is  as  good  a  match  as  physically 
possible.  Sections  of  inverted  magnetization  (slabs)  are  created  by  applying  SLR  inversion 
pulses  in  the  presence  of  slice  selection  gradients.  The  inversion  slab  thickness  is  controlled  by 


Figure  3.  Spin-tagging  sequence.  TR  is  fixed  at  2.7  sec.  In  each  TR  interval,  an  inversion  RF  pulse  is  applied 
first,  followed  by  phase  encoding  steps  starting  at  zero  k  space.  The  time  Thr  is  changed  by  changing  the 
pulse  sequence  control  variable  t_befiv.  The  RF  pulses  are  magnified,  and  t_befiv  and  Tltr  are  not  to  scale 
tactuallv  loneer  than  .shownl. 


12 


changing  the  gradient  strength.  By  applying  the  SLR  pulse  without  an  associated  slice  selection 
gradient,  a  non-selective  slab  containing  inverted  spins  in  the  entire  sensitive  volume  of  the 
transmit  RF  coil  is  created.  By  applying  the  pulse  with  the  gradient,  a  selective  slab  of  specific 
thickness  is  created.  3.  The  so-called  odd-number  interleaved  fc-space  acquisition  [29,30]  has 
been  implemented  to  allow  high  spatial  resolution  in  the  phase  encoding  direction.  With  this 
acquisition,  temporal  resolution  can  be  traded  for  image  spatial  resolution.  Overall,  the  time 
required  to  acquire  an  arterial  spin  tagging  image  at  a  single  slice  is  2  minutes  30  seconds.  This 
longer  time  than  originally  proposed  is  due  to  the  fact  that  images  at  a  multitude  of  TI  times  were 
needed  to  accurately  estimate  the  contribution  of  blood  flow  the  apparent  longitudinal  relaxation 
time,  and  thus  to  estimate  perfusion. 

Our  protocol  using  the  pulse  sequences 

The  protocol  has  been  modified  from  that 
described  in  the  proposal,  to  account  for  the 
longer  scan  time  required  for  the  spin  tagging 
sequence.  For  each  subject,  10  or  more  equally 
spaced  slices  in  the  breasts  are  analyzed. 
Breasts  are  scanned  bilaterally.  A  dual  phased 
array  breast  imaging  RF  coil  (Medical 
Advance,  Inc.,  consigned  specifically  for  this 
project)  is  used  for  the  breast  studies.  Spin 
tagging  measurement  is  based  on  three 
conditions  during  image  acquisition:  1. 
selective  inversion  slab  on  (Selective 
Condition),  2.  non-inversion  slab  on  (Non- 
selective  Condition),  and  3.  no  inversion  slab 
(Regular  Condition).  At  each  slice,  a  time 
series  of  61  images  are  acquired  with  a  fixed 
time  of  repetition  (TR)  of  2.7  seconds,  with  a 
total  scan  time  of  about  3  minutes.  At  each 
time  series.  Conditions  1,  2  and  3  are  applied 
in  sequence.  For  either  the  Selective  or  Non- 
selective  condition,  there  are  seven  TI  cycles 
with  seven  different  TTs  (interval  between  the 
spin  inversion  and  the  effective  center  of  k  space  of  each  image  acquisition).  In  each  TI  cycle, 
four  images  are  acquired.  For  the  Regular  Condition,  five  images  are  acquired  (Figure  4). 
Repeated  image  acquisition  with  the  identical  parameters  is  done  to  improve  the  signal-to-noise 
ratio  (SNR).  The  total  scan  time  for  the  whole  scan  of  12  slice  locations  is  about  30  minutes. 

Contributions  to  data  analysis 

A  data  analysis  technique  described  below  has  been  developed  from  scratch  (no  prior  description 
in  the  literature)  to  properly  interpret  the  timeseries  of  signal  changes  obtained  at  each  pixel.  It 
has  been  developed  taking  into  account  the  characteristics  of  2D  fast  SPGR  data  acquisition.  The 
analysis  rests  on  the  assumption  that  the  2D  fast  SPGR  data  acquisition  does  not  disturb  the 
recovery  curve  of  the  inverted  spins  or  the  equilibrium  states  of  non-inverted  spins.  Although  not 
strictly  true,  the  deviation  from  reality  does  not  significantly  affect  the  outcome  of  the 


Time  series  plot 


Figure.4.  Signal  intensity  timeseries  of  one  pixel 
across  61  images  acquired  using  the  spin  tagging 
sequence.  Selective,  Non-selective  and  Regular 
Tagging  Conditions  are  shown.  For  Selective  and 
Non-selective  Conditions,  images  are  acquired  with 
seven  different  TI  s  in  each,  and  four  images  are 
acquired  at  each  77 .  Images  are  acquired  in  a 
without  pausing  with  TR  =  2.7  seconds. 
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measurements.  Using  images  obtained  under  the  “Regular”  and  the  Selective  or  Non-selective 
Tagging  Conditions,  Tu  or  Tu  is  calculated  at  each  pixel  using  a  semi-log  linear  regression 
(Figure  5), 

J_  ^  In  (5,,,  -5,,,(0)-ln(5,,,  -S,„(j))  (1.8) 

T,,  '  TUi)-Th,(j) 

and 

J_  ^  ln(5.„  -5„.(.))-ln(y„,  -5...0-))  (1.9) 

T,.  ■  r/„(0  -  Ti^m 

where  Sreg  =  the  average  steady-state  signal  intensity  at  a  pixel  under  the  Regular  Condition, 
TIbrii)  =  the  ith  effective  inversion  time  (TI)  (used  for  the  ith  image),  SseKO  =  the  steady-state 
signal  intensity  at  a  pixel  under  the  Selective  Condition  using  the  ith  TI,  and  S„o„(i)  =  the  steady- 
state  signal  intensity  at  a  pixel  under  the  Non-selective  Condition  using  the  ith  TI.  According  to 

the  theory  (reviewed  above),  the  value  off /X  is 
the  difference  between  the  inverse  of  the  Ti 
values  as  expressed  in  Eq.(1.7).  A  derivation  of 
the  above  linear  regression  equations  is  given 
in  the  Appendix. 

As  mentioned  previously,  four  images  are 
acquired  for  each  TI  under  the  Selective 
Condition,  four  images  are  required  for  each  TI 
under  the  Non-selective  Condition,  and  five 
images  are  acquired  under  the  Regular 
Condition.  Generally,  we  have  observed  that 
only  the  last  three  repeated  images  are  at  the 
steady  state  for  their  particular  condition.  The 
first  two  images  the  signal  intensity  is 
approaching  the  steady  state.  Therefore,  only 
these  last  three  are  used  for  the  above  linear 
regression.  A  total  of  21  data  points  are 
available  for  the  calculation  of  T^  or  Tin, 
however,  the  three  low  signal  data  points  near 
the  zero-crossing  line  are  unreliable  and  are 
usually  eliminated  (based  on  an  outlier 
analysis).  Finally,  there  are  18  remaining 
points  for  the  calculation  of  the  Tis.  Because 
the  T\  and//A,  calculations  based  on  this  pulse 
sequence  are  very  sensitive  to  structured  and 
random  noise,  in  particular  the  structured  noise  due  to  motion  (pixel  displacements),  two  noise- 
reduction  techniques  have  been  implemented.  In  the  first  technique,  the  resolution  of  each  image 
is  reduced  from  256  x  256  to  a  lower  resolution,  for  example,  64  x  64,  by  pixel  averaging.  This 
improves  the  signal-to-noise  per  pixel,  and  reduces  the  effect  of  small  tissue  motions  occurring 
during  the  61  images  being  analyzed.  This  reduction  is  due  to  the  simple  fact  that  any  particular 


Steady-state  signal  intensity 
difference  as  a  function  of  77 


Thr  =  Start  of  image  acquisition,  measured  from  time  magnetization  is 
inverted;  Sreg  =  Average  signal  from  last  three  images  during  Regular 
Tagging  Condition;  S„j  =  Signal  from  steady-state  image  during 
Selective  Tagging  Condition;  Sno„  =  Signal  from  steady-state  image 
during  Non-selective  Tagging  Condition. 

Figure  5.  Example  semi-log  linear  regression  method  for 
calculating  T\  values  (Ju  for  Selective  Tagging 
Condition,  and  Tin  for  Non-selective  Condition)  and 
effective  perfusion  //A,.  The  calculation  is  based  on  the 
number  of  points  remaining  after  use  of  minimum-mean- 
square-error  technique  to  remove  outliers. 
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motion  results  in  a  lower  fractional  change  in  the  pixel  volume.  The  second  technique  uses  a  data 
elimination  strategy  to  reduce  the  mean-square-error  of  the  semi-log  regression  used  for 
determination  of  the  Tis.  This  technique  permits  the  elimination  of  up  to  five  noisy  data  points 
(out  of  18)  to  reduce  the  Ti  regression  error. 

Contributions  to  image  processing 

A  computer  program,  called  BrView,  has  been  developed  for  display,  analysis,  and  interpretation 
of  the  images  obtained  in  the  breast  imaging  protocol,  starting  from  previous  programs 
developed  by  the  PI  for  neuro-functional  MRI  research  (as  described  in  the  proposal).  The 
program  is  written  in  C,  the  X  Window  System,  and  Motif.  For  all  of  the  analysis  methods 
described  above,  the  values  of  Tu,  Tu  and//A,  at  different  pixel  locations  on  each  slice  are 
calculated.  Images,  which  we  call  feature  images,  are  formed  based  on  these  values.  The  three 
most  useful  feature  images  are  those  showing  Tin  values,  the  perfusion  image  showing//A, 
values,  and  the//X,  error  image  based  on  the  standard  deviations  of//A,  values.  These  feature 
images  contain  pixels  with  a  wide  range  of  values.  However,  values  of  T\,f/X  andf/X  errors  of 
normal  breast  tissues  are  within  narrow  ranges,  while  pixels  displacing  large//X,  Ti  above  500 
ms,  but  small//X,  error  are  regarded  as  suspicious  of  malignant  disease.  Pixels  containing  fat 
display  very  short  Ti  and  the  perfusion  values  at  these  pixels  tend  to  be  inconsistent.  To  display 
these  feature  images,  an  appropriate  positive  shift,  scaling  and  window  size  is  applied  [31]. 

This  BrView  program  is  used  to  display  and  interpret  the  feature  images  for  identification  of 
suspicious  regions.  Since  a  sizable  breast  lesion  covers  several  pixels  within  the  same  slice  and 
possibly  neighboring  slices,  these  pixels  should  have  similar  tissue  characteristics.  It  is  assumed 
that  pixels  with  a  breast  lesion  will  have  approximately  equal  Ti  values,  and  that  these  Ti  values 
will  be  relatively  high  compared  to  fat,  and  also  show  high  perfusion  rates  with  low  perfusion 
standard  errors.  When  a  pixel  is  found  to  have  a  moderately  long  Ti,  and  a  high//A,,  and  the//X, 
error  is  low,  surrounding  pixels  within  the  same  slice  and  across  slices  are  also  checked  to  see  if 
they  have  similar  characteristics.  The  “suspicion  index”  for  this  pixel  is  defined  as  the  percentage 
of  surrounding  pixels  also  classified  as  suspicious.  By  applying  this  method  to  all  pixels  in  all 
slices,  “suspicion  index  maps”  are  created.  If  the  pixel  has  a  suspicion  index  exceeding  a 
threshold,  this  pixel  would  be  highlighted  on  the  display.  Clearly,  the  level  of  suspicion  here  is 
qualitative.  Neither  a  rigorous  statistical  analysis  of  the  procedure,  nor  a  statistical  comparison 
and  analysis  of  suspicion  regions  with  the  results  of  dynamic  contrast  enhanced  MRI,  has  not 
been  carried  out  yet.  Nevertheless,  the  goal  is  to  develop  a  quantitative  tool  for  breast  cancer 
detection.  The  suspicion  index  threshold  for  identification  of  malignant  disease,  as  well  as 
independent  thresholds  for/ZX,  T\  and  the  standard  deviation  of f/X,  should  ultimately  be  based 
on  the  results  of  a  large  set  of  patient  studies. 

Based  on  the  small  set  of  studies  that  have  been  done,  the  thresholds  used  to  identify  suspicious 
pixels  are  as  following: 

Ti:  >  0.5  sec 

f/X:  >0.1  sec'* 

STD  of  f/X:  <  0.1  sec'*. 

and  the  suspicion  index  threshold,  used  for  asserting  a  “high”  probability  that  the  tissue  is 
abnormal,  is  set  at  20  %.  A  quantitative  analysis  of  the  false  positive  and  negative  error  rates 
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associated  with  this  threshold  has  not  been  done.  The  BrView  program  allows  clinician 
researchers  to  change  the  thresholds  and  observe  the  changes  in  the  suspicion  index  maps.  This 
visualization  program  can  simultaneously  display  feature  images  and  standard  clinical  MR 
images.  Parameters  T\,flX,  and  the  standard  deviation  of //A,  at  any  pixel  location  of  the  breast 
objects  can  be  found  by  a  click  of  the  mouse  on  that  location.  If  a  dynamic  contrast  enhancement 
study  has  been  done,  a  click  of  the  mouse  reveals  a  plot  of  the  change  of  signal  intensity  as  a 
function  of  time  at  this  location.  The  suspicious  regions  can  also  be  mapped  to  high-resolution 
traditional  clinical  images  (256  x  256),  which  are  always  acquired  for  comparison  during  each 
subject  study.  This  mapping  provides  a  radiologist  the  traditional  view  of  the  breast  images  with 
the  suspicious  regions  identified.  In  summary,  BrView  is  designed  to  allow  pixel  cross¬ 
registration  and  comparison  between  results  from  different  MR  imaging  techniques.  It  provides 
an  efficient  working  environment  for  clinicians  making  diagnostic  decisions. 


Other  image  processing  contributions  are  methods  to  reduce  tissue  misregistration  that  cause 
inaccuracy  in  the  T\  estimation.  Breast  images  are  expected  to  precisely  align  in  all  61  images 


Reference  image 

Image  to  be  realigned 

Region  for  local  cross  correlation 
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Figure  6.  Local  cross  correlation  is  performed  Figure  7.  Kidney  dialysis  phantom  flow  study:  Arterial  spin- 
between  the  binary  image  and  reference  binary  tagging  was  applied  at  four  different  flow  rates:  0,  0.0222, 
image.  The  maximum  correlation  occurs  when  the  0.0640  and  0.0787  ml/sec/cm^.  The  corresponding  calculated// 
two  objects  maximally  overlap.  X  values  were  0.0251,  0.138,  0.331  and  0.409  sec  ’. 

acquired  in  the  timeseries  at  each  slice  location.  Unfortunately,  the  patient  may  move  during 
the  two  minute  30  second  scan,  resulting  in  a  change  in  breast  position  and  orientation.  The 
motion  estimation  used  here  is  based  on  binary  segmentation,  following  by  statistical  best-fit 
alignment  of  defined  regions  across  the  collection  of  61  images.  Geometric  constraints  of 
patient  motion  and  breasts,  within  the  scanner  and  dedicated  RF  coil,  suggest  that  only  lateral 
displacements  in  the  X  and  Y  directions  should  be  considered.  The  studies  performed  suggest 
that  typical  displacements  are  not  large,  within  1  cm.  Such  motions  are  correctable  with 
segmentation  and  region  alignment  (see  Figure  6). 


Results  of  in-vitro  and  in-vivo  testing 

The  success  of  the  arterial  spin  tagging  pulse  sequence  in  measuring  perfusion  was  first  proven 
by  the  result  of  flow  phantom  studies.  All  studies  showed  that  there  was  a  linear  relationship 
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between//A,  and  actual  flow.  All  evaluated  flows  were  within  the  range  of  tissue  perfusion 
(Figure  7).  Furthermore,  the  Tin  values  obtained  in  the  non-selective  condition,  at  four  different 
flow  rates,  were  as  expected,  specifically  1.232  +/-  0.034  sec.  A  brain  imaging  study  was 
performed  to  further  evaluate  the  arterial  spin  tagging  pulse  sequence.  In  1 1  different  regions  of 
interest  (ROIs),  defined  at  several  different  slice  locations,  white  matter  had  T\n  equal  to  0.495 
+/-  0.021  sec,  Jis  equal  to  0.487  +/-  0.022  sec,  and//X  equal  to  0.030  +/-  0.075  sec  ranging 
from  0  to  0.131  sec  Gray  matter  in  this  study  had  Tin  equal  to  0.665  +/-  0.065  sec,  Tis  equal  to 
0.629  +/-  0.035  sec,  and//A,  equal  to  0.081  +/-  0.054  sec  ranging  from  0.037  to  0.170  sec  In 
two  normal  breast  studies,  at  14  different  ROIs  defined  at  several  different  locations, 
parenchyma  (glandular  tissue)  had  Tin  equal  to  0.772  +/-  0.088  sec,  T^  equal  to  0.709  +/-  0.062 
sec,  and//X  equal  to  0.109  +/-  0.053  sec  ranging  from  0.043  to  0.190  sec’*.  Normal  fat  tissue 
had  Tin  equal  to  0.377  +/-  0.069  sec,  T^  equal  to  0.363  +/-  0.056  sec,  and  fiX  equal  to  0.086  +/- 
0.073  sec  ■*,  ranging  from  0.003  to  0. 177  sec  The  values  resulting  from  the  head  and  normal 
breast  studies  were  within  expected  ranges. 

Results  of  patient  studies 

Introduction 

For  this  and  all  cases  summarized  below,  suspicious  pixels  were  identified  based  on  the 
following  criteria: 

Ti:  >  0.5  sec 
//X,:  >0.1  sec’* 

STD  of//X:  <0.1  sec’*. 

and  the  suspicion  index  threshold  was  set  at  20.2  %.  Analysis  was  done  using  the  program 
BrView. 


Case  1 

Case  Study  #1:  palpable  mass 
Date  of  study:  Dec  10,  1997 

Level  of  suspicion  based  on  Tin  and  perfusion:  Moderately  High 
Clinical  MR  images:  Negative  (lesion  not  visualized) 

Assessment  based  on  spin  tagging:  Positive  for  malignancy 
Biopsy  result:  benign  (spin  tagging  INCONSISTENT  with  biopsy) 

Discussion:  Figure  8  shows  the  three  feature  images,  and  an  image  acquired  under  the  “Regular 
Tagging  Condition”  of  the  arterial  spin  tagging  sequence.  All  three  feature  images  are  64  x  64. 
The  Ti  image  has  excellent  SNR  and  clearly  shows  the  distribution  of  Ti  values.  In  particular, 
bright  areas  at  the  centers  of  the  breasts  represent  glandular  tissue  with  long  Ti  values.  The 
corresponding  perfusion  image  has  relatively  low  SNR.  Some  pixels  at  the  edge  of  the  breasts 
have  high  intensity  due  to  high  calculation  error,  not  high//X.  The  calculation  at  edge  pixels  is 
extremely  sensitive  to  breast  motion,  as  well  as  to  the  presence  of  Gibbs  ringing  artifacts 
common  in  low  resolution  images.  Pixels  with  high  T\  values  generally  have  less  perfusion  error. 
Because  the  perfusion  estimates  are  more  sensitive  to  errors  in  Ti  when  Ti  is  low.  Fortunately, 
this  property  of  the  error  should  not  limit  the  identification  of  malignant  tissue  very  much, 
because  malignant  tissue  will  typically  have  relatively  high  Ti  values  (>  500  ms).  Figure  9 
shows  the  overlay  of  the  suspicious  regions  onto  the  high-resolution  clinical  images.  Highlighted 
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(pure  white)  pixels  on  the  left  breast  (right  side  of  image)  on  slices  4,  5, 6  and  7  may  be 
identifying  one  contiguous  mass.  Physical  examination  of  this  patient  revealed  a  mass  at  the  nine 
o'clock  position  in  the  left  breast,  consistent  with  the  suspicious  cancer  region  identified  by  spin 
tagging. 

Case  2 

Case  Study  #2:  palpable  mass 

Date  of  study:  Nov  16,  1997 

Level  of  suspicion  based  on  Tin  and  perfusion:  Low 

Clinical  MR  images:  Positive  (lesion  visualized,  consistent  with  palpation) 

Assessment  based  on  spin  tagging:  Negative  for  malignancy 
Biopsy  result:  benign  (spin  tagging  CONSISTENT  with  biopsy) 

Discussion:  The  low  suspicion  index,  based  on  T\  relaxation  time  and  perfusion,  at  and  around 
the  palpated  lesion  location  was  strong  indication  that  the  lesion  was  benign.  Only  one 
suspicious  pixel  was  found  at  that  location  at  the  20.2%  threshold.  It  was  judged  that  the  breasts 
did  not  contain  malignant  tissue.  This  case  is  shown  in  Figures  10(a)-(d). 

Case  3 

Case  Study  #3:  Palpable  mass 

Date  of  study:  Aug  6,  1997 

Level  of  suspicion  based  on  Tin  and  perfusion:  High 

Clinical  MR  images:  Negative  (mass  not  clearly  visualized,  prior  cyst  visualized) 

Result  based  on  spin  tagging:  Positive  for  malignancy 
Biopsy  result:  Positive  (spin  tagging  CONSISTENT  with  biopsy) 

Discussion:  The  high  suspicion  index,  based  on  T\  relaxation  time  and  perfusion,  at  and  around 
the  lesion  detected  by  palpation,  was  strong  indication  that  the  lesion  was  malignant.  This  case  is 
shown  in  Figures  ll(a)-(d). 

Case  4 

Case  Study  #4:  Palpable  mass 

Date  of  study:  Aug  7,  1998 

Level  of  suspicion  based  on  Tin  and  perfusion:  High 

Clinical  MR  images:  Positive  (lesion  visualized) 

Assessment  based  on  spin  tagging:  Positive  for  malignancy 

Assessment  based  on  dynamic  contrast-enhanced  study:  Positive  for  malignancy 

Biopsy  result:  Negative  (spin  tagging  NOT  CONSISTENT  with  biopsy) 

Discussion:  High  suspicion  indices,  based  on  relaxation  time  and  perfusion,  were  found  in 
several  breast  regions.  The  visible  lesion  location  on  clinical  MRI  was  consistent  with  that  of  the 
palpable  mass.  It  overlapped  but  was  considerably  smaller  than  that  found  by  the  suspicion 
indices.  First-pass  contrast  enhancement  study  results  correlated  strongly  with  the  high  perfusion 
regions  detected  by  spin  tagging.  Nevertheless,  biopsy  showed  the  mass  to  be  benign.  This  case 
is  shown  in  Figures  12(a)-(d). 

In  summary,  there  was  only  partial  agreement  (2  of  4)  between  spin  tagging  and  biopsy 
assessments  of  malignant  disease.  More  patient  studies  (as  planned  in  project),  as  well  as  follow¬ 
up  of  patients,  is  needed.  Spin  tagging  may  be  detecting  abnormal  regions  not  evaluated  by  the 
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biopsies.  Arterial  spin  tagging  results  were,  subjectively,  more  convincing  if  the  suspicious 
regions  overlapped  abnormal  regions  seen  on  standard  clinical  T\,T2  and  proton-density 
weighted  images.  Unfortunately,  standard  clinical  images  often  did  not  reveal  any  abnormal 
tissue,  concurring  with  the  view  that  standard  clinical  MRI  is  not  definitive  in  the  assessment  of 
breast  lesions.  Dynamic  first-pass  contrast-enhanced  imaging  is  definitely  needed  to  help 
validate  the  spin  tagging  technique,  and  this  imaging  will  be  done  on  all  subsequent  patients  in 
this  project.  In  addition,  results  from  X-ray  mammography  and  ultrasonography  will  be 
compared,  as  has  been  done  for  prior  studies  (see  page  40). 

Future  Implementations 

Arterial  spin  tagging  with  an  echo  planar  type  data  acquisition 

To  derive  the  semi-log  regression  analysis,  it  was  assumed  that  SPGR  data  acquisition  would  not 
disrupt  the  recovery  of  longitudinal  magnetization  of  inverted  spins,  nor  the  equilibrium  of  the 
non-inverted  spins.  Unfortunately,  this  assumption  is  a  cause  of  error  in  the  method,  and  our 
simulations  have  shown  that  the  T\  calculated  based  on  the  assumption  was  slightly  lower  than 
the  true  T\.  Using  the  simulation  results,  a  conversion  table  was  developed  to  convert  the 
estimate  derived  based  on  the  assumption  to  a  better  T1  estimate.  Nevertheless,  a  better  approach 
would  be  to  replace  the  2D  fast  SPGR  data  acquisition  scheme  with  an  echo  planar  data 
acquisition  scheme.  Echo  planar  imaging  (EPI),  requires  only  a  single  RF  pulse  for  acquisition  of 
the  entire  data  set  for  one  image.  The  disturbance  due  to  the  repeated  RF  pulses  to  the  recovery 
curves  of  the  inverted  spins  or  the  equilibrium  states  of  the  non-inverted  spins  would  be  much 
more  predictable.  Also,  EPI  is  faster  than  2D  fast  SPGR,  so  the  77  cycles  could  be  redesigned, 
and  a  multi-slice  technique  introduced,  to  shorten  the  scan  time  by  a  factor  of  4-8.  This 
implementation  would  make  possible  high  resolution  scanning  of  the  entire  breast  tissue  in  under 
10  minutes.  Starting  September  1998,  the  UC  Davis  Medical  Center  will  operate  a  GE  Signa 
Horizon  Echospeed  MR  system.  This  system  will  be  capable  of  echo  planar  imaging  (128  x  128 
single  shot  is  typical,  but  up  to  256  x  256  single  shot  can  be  done).  Conversion  of  the  arterial 
spin  tagging  sequence  to  incorporate  echo  planar  acquisition  will  be  carried  out.  In  addition, 
conversion  will  also  be  done  for  the  existing  spin  tagging  sequence  that  is  based  on  2D  fast 
SPGR. 

Magnetization  transfer  implementation 

The  second  major  assumption,  used  in  the  derivation  of  the  perfusion  formula,  is  that  the 
contribution  of  the  magnetization  transfer  term  in  Eq.  (1.1)  is  negligible.  To  allow  inclusion  of 
this  term,  the  spin  tagging  pulse  sequence  must  include  a  binomial  or  high  order  RF  pulse  to 
prepare  the  spins  prior  to  data  acquisition.  Differential  magnetization  transfer  effects  have  been 
detected  in  head  and  neck  neoplasms  [32].  Development  of  this  pulse  sequence  modification  to 
estimate  the  size  of  the  magnetization  transfer  contribution  will  be  completed. 
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Figure  9.  Dec  10,  1997  study.  The  suspicious  pixels  (threshold  20.2%),  showing  moderate  Tj  ,  h\ghf/X,  \ovj  f/X  error,  are 
identified  on  six  (a-f)  contiguous  slices,  displayed  on  fat  suppressed  weighted  fast  spin  echo  images. 
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Figure  10(d).  Nov.  16,  1997  study.  The  T,  and  f/X  values  of  pixels  representing  the 
lesion  identified  by  proton  density  and  T2  images,  are  superimposed  on  the  histogram 
of  the  Tj  andf/X  values  from  the  entire  breast.  Shows  that  criteria  for  suspicuous 
pixels  is  based  on  moderate  Tj  and  high/7A.  See  text  for  complete  criteria. 
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Figure  11(b).  Aug  6,  1997  study.  Feature  images  (top  row),  suspicion  index  map  (bottom  left,  threshold  20.2%) 
and  mapping  (pure  white  pixels)  on  proton  density  and  T2  images,  (bottom  row).  Also  shown:  Axial  T^. 
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density  and  T2  images.  Also  shown:  Axi 
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Figure  11(d).  Aug  6,  1997  study.  The  T,  and  f/X  values  of  pixels  representing  the  lesion  identified 
by  proton  density  and  T2  images  in  Figure  1 1(c),  are  superimposed  on  the  histogram  of  the  Tj  and^A 
values  from  the  entire  breast.  Shows  that  criteria  for  suspicion  is  based  on  moderate  and  high //A. 
See  text  for  complete  criteria. 
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i  study.  Feature  images  (top  row),  suspicion  index  map  (bottom  left,  threshold  20.2^ 
images  are  shown  without  suspicion  index  map  (bottom  row).  Also  shown:  Axial  J, 
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(7)  CONCLUSIONS 

In  this  year  of  the  project,  we  have  developed  the  necessary  pulse  sequences  for  anatomical  and 
functional  imaging  of  breast  tissue,  and  the  necessary  software  for  the  analysis,  display  and 
interpretation  of  this  data.  These  developments  will  allow  us  to  evaluate  complete  breast  imaging 
studies  that  are  composed  of  anatomical  scans,  dynamic  first  pass  contrast  enhanced  scans,  and 
arterial  spin  tagging  scans.  The  integrated  nature  of  the  software  allows  easy  and  consistent 
interpretation  of  both  spatial  dependencies  as  well  as  temporal  dependencies  of  signal  changes 
indicative  of  disease.  The  in-vitro  and  in-vivo  studies  done  thus  far  indicate  that  first  and 
foremost,  arterial  spin  tagging  is  a  sensitive  and  reliable  method  for  measuring  T\  and  parameters 
related  to  tissue  perfusion.  Second,  they  indicate  that  arterial  spin  tagging  is  a  viable  alternative 
to  dynamic  first-pass  contrast-enhanced  imaging.  The  fact  that  arterial  spin  tagging  can  be  easily 
added  to  any  standard  breast  imaging  protocol,  without  requiring  special  nursing  or  MR 
technologist  expertise,  means  that  it  is  especially  attractive  for  routine  clinical  use.  It  could  be 
used  to  supplement  contrast-enhanced  studies  in  the  initial  study  of  the  patient,  and  used 
exclusively  in  follow-up  studies.  Since  10/1/98,  six  subjects  have  been  studied  without  the 
dynamic  first-pass  contrast-enhanced  scan,  and  two  have  been  studied  with  this  scan. 

Comparison  of  contrast-enhanced  images  versus  arterial  spin-tagged  images  has  not  yet  been 
made.  In  particular,  no  statistical  comparison  of  the  dynamic  study  results  (considered  to  be  the 
gold  standard  for  non-invasive  assessment  of  malignant  disease)  and  the  arterial  spin  tagging 
results  have  been  performed.  However,  there  is  considerable  spatial  overlap  of  suspicious  regions 
obtained  by  these  methods  noted  on  visual  inspection.  Collecting  data  on  human  subjects,  and 
statistical  confirmation  of  this  agreement,  will  be  major  focuses  of  next  year’s  work. 
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(9)  APPENDIX 


Contributions  to  Data  Analysis 

The  regression  method  derived  in  this  project  resolves  a  key  confound  that  is  present  in  all 
arterial  spin  tagging  sequences,  that  of  the  mismatch  of  the  slice  profiles  of  the  inversion  and 
excitation  RF  pulses.  Our  sequence  is  longer  than  many  arterial  spin  sequences  reported  in  the 
literature,  but  we  are  quite  sure  that  these  other  methods  fail  to  account  for  the  mismatch,  and 
result  in  dubious  measurements  of  perfusion.  It  is  interesting  to  note  that  the  mismatch  results  in 
signal  contamination  who  spatial  dependence  resembles  an  inverted  Ti  weighted  image,  thus, 
they  are  easily  mistaken  for  true  perfusion  maps  (i.e.  in  the  brain,  CSF  is  bright,  followed  by 
gray  matter,  then  white  matter,  just  as  true  perfusion  would  present). 

Derivation  of  semi-log  linear  regression  method 

In  this  derivation,  it  is  assumed  that,  because  the  excitation  RF  flip  angle  is  small,  repeated 
application  of  these  pulses  for  data  acquisition  does  not  disturb  Ti  recovery  of  the  inverted  spins, 
nor  the  equilibrium  magnetization  of  the  imaged  spins  if  not  inverted.  Given  this  assumption, 
under  a  spin  tagging  condition,  the  steady-state  longitudinal  magnetization  at  the  end  of  each  TR 
period  just  before  the  next  inversion  pulse  would  be 


where,  Mtoss  =  the  steady-state  longitudinal  magnetization  at  the  end  of  each  TR  period.  Mo  =  the 
longitudinal  magnetization  under  fully  relax  condition,  TR  =  the  repetition  time,  and 


f  1 


for  the  Selective  Tagging  Condition 


for  the  Non-selective  Tagging  Condition 


Therefore, 


^tOss 


=  Mn 


\+e~‘’  ™ 


(A.1) 


The  longitudinal  magnetization  at  time  t  after  the  spin  inversion  at  steady-state  is 
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(A.2) 


=  Mo(1-7-^) 

1  +  e 

As  a  result,  the  signal  measured  at  echo  time  TE  is 


5,=Mo(l-— ^)sin(0> 
l  +  e 


(A.3) 


where,  77  =  the  time  of  inversion  from  the  inversion  pulse  to  the  effective  center  of  k  space,  TE  ■■ 
the  time  of  echo  at  each  read-out  period,  and  6  =  the  excitation  flip  angle. 

Under  the  Regular  Tagging  Condition,  the  longitudinal  magnetization  Mreg{t)  at  any  time  t  stays 
at  equilibrium.  Specifically, 


The  signal  measured  at  echo  time  TE  becomes 


5,,  =MoSin(6/)e 


(A.4) 


(A.5) 


Thus, 


S,  2e~"” 

^  l  +  e'*™ 


(A.6) 


Thus,  considering  the  ith  inversion  time  (TI): 

5,,, -5,(0 

S„-S,Q)  e'*""’ 


_  g-6(n(/)-r/(/)) 


(A.7) 


where,  7/(0  =  the  time  of  inversion  at  the  ith  TI  cycle,  5,(0  =  the  signal  at  a  pixel  when  TI{i)  is 
used,  7/(0  =  the  time  of  inversion  at  the  jtb  TI  cycle,  5,(/)  =  the  signal  at  a  pixel  when  TI(j)  is 
used.  Taking  the  logarithm  of  both  sides  yields. 


(  I 


In  -  5,  (/))  -  In  -  5,  (/))  =  -&(77(/)  -  TI  (/)) 

and  solving  for  b  gives, 

^  \n{S,,^-S,(i))-\HS,,^-S,{j))  (A.8) 

Tm-TI(j) 

The  above  equation  shows  the  semi-log  linear  regression  relationship  between  (Sreg  -  St)  and  TI. 
The  formula  can  be  applied  to  both  Selective  and  Non-selective  tagging  conditions,  do  derive  bs 
and  b„,  respectively.  The  exact  values  of  the  effective  TTs  are  uncertain  (due  data  collection  is 
spread  out  over  1  second).  Nevertheless,  it  is  reasonable  to  assume  that  these  values  each  have 
the  same  amount  of  time  error  relative  to  their  true  values.  Because  TI(i)  and  TI(j)  is  subtracted 
in  the  denominator,  the  uncertainty  of  the  timing  of  this  “effective”  TI  cancels  out. 

Elimination  of  slice  profile  mismatch  effects 

The  semi-log  relationship  eliminates  the  errors  due  to  mismatch  between  the  slice  profiles  of  the 
inversion  and  data  acquisition  RF  pulses.  Because  the  RF  pulses  used  for  inversion  and  data 
acquisition  are  different  under  the  Selective  Condition,  (The  Non-selective  Condition  does  not 
have  this  problem  because  the  RF  inversion  pulse  inverts  spins  everywhere.). 

The  geometric  mismatch  between  the  inversion  and  read-out  profiles  would  lead  to  the 
contamination  of  signal  from  non-inverted  spins.  Thus  the  apparent  signal  from  the  Selective 
Condition  can  be  modeled  as 

=  (1  -  +  err  (A.9) 

where  Ssei  -  the  signal  of  a  pixel  if  no  slice  profile  mismatch,  =  the  signal  of  a  pixel  actually 
measured,  and  err  =  the  amount  of  percent  error  due  to  slice  profile  mismatch.  Then, 

-^  =  {l-err)^+err 
S  s 

^  rfio  ^  rjta 


and 


^sel„r 


Ssel  _  ^reg 


-err 


Sreg  ^-err 


Thus, 
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l-err- 


(A.10) 


+  err 


\-err 


\-err 


Comparing  with  Eq.  (A.7),  we  have 


^reg  (0  _  (TI(iyTI(i)) 

^reg  -  (/) 


where, 


b. 


1  f 

— h  —  =  the value  under  the  Selective  Condition 

7;  X 


Thus, 


*’■ - ¥I(iPFm - 

Eq.  (A.  1 1)  corresponds  exactly  as  Eq.  (A.8).  In  other  words,  by  using  the  semi-log  linear 
regression  method,  the  problem  due  to  geometric  mismatch  between  the  inversion  and  read-out 
slice  profiles  is  eliminated. 

Phantom,  normal  subject,  and  patient  study  notes 

These  notes  (pages  41-50,  presented  in  two-column  landscape  mode)  were  written  immediately 
after  each  study,  and  provide  brief  subject  history,  and  detailed  listing  of  pulse  sequences  and 
parameters  used  in  the  protocol,  and  problems  noted.  Preliminary  accounts  of  studies  prior  to 
1 1/19/97  have  been  reported  [33,34,35]. 


Inventory  of  Breast  Imaging  Studies  to  Date 


Test  Phantoms 

Normal 

Subject 

Patient 
(no  CE-MRI) 
No  biopsy 

Patient 
(no  CE-MRI) 
Biopsy 

Patient 

(with  CE-MRI) 
Biopsy 

Prior  to  10/1/98 

9 

4 

3 

3b,  2m 

0 

Since  10/1/98 

5 

1 

1 

4b 

2b 

CE-MRI  (dynamic  first-pass  contrast-enhanced  MRI).  Number  preceding  b  indicates  number  of  subjects  in  which 
the  biopsy  revealed  benign  disease.  Number  preceding  m  indicates  number  of  subjects  in  which  the  biopsy  revealed 
malignant  disease.  Patients  without  biopsy  did  not  meet  clinical  criteria  and  mass  was  assumed  benign. 
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Date  of  study:  5/11/1997 
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001:  £gtl^ew3b,  default  setting 

of  study:  7/6/97  autoprescan:  Rl  =  6,  R2  =15,  TG  =  109 
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3D  scan  coronal  cut  (Series  4):  series  7:  location:  a38,  R1  =  6  ,  R2  =  15  ,  TG  =  106 


Place  the  Meat  phantom  with  the  7-up  phantom 
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Number  of  images  =  10,  br_norm  =  2,  br_tino  =  2,  br_acqs 
br_TI  =  100  ms,  ti_inc  =  200  ms,  br_TR  =  2.7  s 
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locates  at  left  breast,  6  o'clock,  about  1/3  in  from  the  nipple.  2.  2D  Past  Spin  Echo  (Series  2)  FOV 
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to  have  a  better  fit  of  the  foam  cast,  tighter  strap, 

love,  and  to  breath  smoothly  with  stomach.  4.  3D  scan  (scan  #2  in  protocol,  2‘37")  with  flip  angle  20  degree,  (series  4) 

28  slices  with  the  center  at  A33.6. 
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REPLY  TO 
ATTENTION  OF 


DEPARTMENT  OF  THE  ARMY 
US  ARMY  MEDICAL  RESEARCH  AND  MATERIEL  COMMAND 
504  SCOTT  STREET 
FORT  DETRICK,  MD  21702-5012 


MCMR-RMI-S  (70-ly) 


1  July  03 


MEMORANDUM  FOR  Administrator,  Defense  Technical  Information 
Center  (DTIC-OCA) ,  8725  John  J.  Kingman  Road,  Fort  Belvoir, 
VA  22060-6218 

SUBJECT:  Request  Change  in  Distribution  Statement 


1.  The  U.S.  Army  Medical  Research  and  Materiel  Command  has 
reexamined  the  need  for  the  limitation  assigned  to  techriical 
reports  written  for  this  Command.  Request  the  limited 
distribution  statement  for  the  enclosed  accession  numbers  be 
changed  to  "Approved  for  public  release;  distribution 
unlimited. "  These  reports  should  be  released  to  the  National 
Technical  Information  Service. 


2.  Point  of  contact  for  this  request  is  Ms.  Kristin  Morrow  at 
DSN  343-7327  or  by  e-mail  at  Kristin.Morrow@det.amedd.army.mil 


FOR  THE  COMMANDER; 


RINEfi^T 
ief  of  Staff  for 
'at ion  Management 


End 
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